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ABSTRACT 

Radiation damage to space-based Charge-Coupled Device (CCD) detectors creates 
defects which result in an increasing Charge Transfer Inefficiency (CTI) that causes 
spurious image trailing. Most of the trailing can be corrected during post-processing, 
by modelling the charge trapping and moving electrons back to where they belong. 
However, such correction is not perfect - and damage is continuing to accumulate in 
orbit. To aid future development, we quantify the limitations of current approaches, 
and determine where imperfect knowledge of model parameters most degrade mea¬ 
surements of photometry and morphology. 

As a concrete application, we simulate 1.5 x 10® “worst case” galaxy and 1.5 x 10® 
star images to test the performance of the Euclid visual instrument detectors. There 
are two separable challenges: If the model used to correct CTI is perfectly the same as 
that used to add CTI, 99.68 % of spurious ellipticity is corrected in our setup. This is 
because readout noise is not subject to CTI, but gets over-corrected during correction. 
Second, if we assume the first issue to be solved, knowledge of the charge trap density 
within Ap/p = (0.0272 ±0.0005)%, and the characteristic release time of the dominant 
species to be known within Ar/r = (0.0400 ± 0.0004)% will be required. This work 
presents the next level of definition of in-orbit CTI calibration procedures for Euclid. 

Key words: space vehicles: instruments — instrumentation: detectors — methods: data analysis 


1 INTRODUCTION 

The harsh radiation environment above the Earth’s atmo¬ 
sphere gradually degrades all electronic equipment, includ¬ 
ing the sensitive Charge-Coupled Device (CCD) imaging de¬ 
tectors used in the Hubble Space Telescope (HST) and Gaia 
2008[ ), and proposed for use by Euclid | |Lau-| 
). CCD detectors work by collecting pho¬ 
toelectrons which are stored within a pixel created by an 
electrostatic potential well. After each exposure these elec¬ 
trons are transferred via a process called clocking, where 
alternate electrodes are held high and low to move charge 
through the pixels towards the serial register. The serial 
register is then clocked towards the output circuit where 


charge-to-voltage conversion occurs providing an output sig¬ 
nal dependent on the charge contained within a pixel. The 
amount of charge lost with each transfer is described by 
the Charge Transfer Inefficiency (CTf). One of the results 
of radiation-induced defects within the silicon lattice is the 
creation of charge traps at different energy levels within the 
silicon band-gap. These traps can temporarily capture elec¬ 
trons and release them after a characteristic delay, increas¬ 
ing the CTf. Any electrons captured during charge transfer 
can re-join a charge packet later, as spurious charge, often 
observed as a charge tail behind each source. 

Charge trailing can be (partially) removed during im¬ 
age postprocessing. Since charge transfer is the last process 


([Lindegren et al. 


reijs et al. 2011 
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to happen during data acquisition, the fastest and most suc¬ 
cessful attempts to correct CTI take place as the second step 
of data reduction, right after the analogue-digital converter 
bias has been subtracted, {e.g. Bristow|2003 1. By modelling 
the solid-state physics of the readout process in HST’s Ad¬ 
vanced Camera for Surveys (ACS), then iteratively reversing 
the model, Massey et al. ( |2010| demonstrated a 10-fold re¬ 
duction in the level of charge trailing. The algorithm was 
sped up by Anderson & Bedin (20101 and incorporated into 
STScI’s HST default analysis pipeline (Smith et al. 20121 . As 
the radiation damage accumulated, the trailing got bigger 
and easier to measure. With an updated and more accurate 
HST model, Massey (20101 achieved a 20-fold reduction. In 
an independent programme for Gaia, Short et al. (20131 
developed a model using different underlying assumptions 
about the solid-state physics in CCDs. Massey et al. (20141 
created a meta-algorithm that could reproduce either ap¬ 
proach through a choice of parameters, and optimised these 
parameters for HST to correct 98% of the charge trailing. 

The current level of achievable correction is acceptable 
for most immediate applications. However, radiation dam¬ 
age is constantly accumulating in HST and Gaia-, and in¬ 
creasing accuracy is required as datasets grow, and statisti¬ 
cal uncertainties shrink. One particularly challenging exam¬ 
ple of stringent requirements in future surveys will be the 
measurement of faint galaxy shapes by Euclid. 

In this paper, we investigate the effect of imperfect CTI 
correction, on artificial images with known properties. We 
add charge trailing to simulated data using a CTI model 
M, then correct the data using a CTI model with imper¬ 
fectly known parameters, M -|- 5M. After each stage, we 
compare the measured photometry (flux) and morphology 
(size and shape) of astronomical sources to their true (or 
perfectly-corrected) values. We develop a general model to 
predict these errors based on the errors in CTI model pa¬ 
rameters. We focus on the the most important parameters 
of a ‘volume-driven’ CTI model: the density pi of charge 
traps, the characteristic time Ti in which they release cap¬ 
tured electrons, and the power law index j3 describing how 
an electron cloud hlls up the physical pixel volume. 

This paper is organised as follows. In Sect. we sim¬ 
ulate Euclid images and present our analysis methods. In 
Sect. we address the challenge of measuring an average 
ellipticity in the presence of strong noise. We present our 
CTI model and measure the CTI effects as a function of trap 
release timescale r in Sect.|^ Based on laboratory measure¬ 
ments of an irradiated CCD273 (Endicott et al. 20121, we 
adopt a baseline trap model M for the Euclid VIS instru¬ 
ment (Sect.|^. In this context, we discuss how well charge 
trailing can be removed in the presence of readout noise. 
We go on to present our results for the modified correction 
model (M -I- dM) and derive tolerances in terms of the trap 
parameters based on Euclid requirements. We discuss these 
results in Sect. [ 6 ] and conclude in Sect.[3 


2 SIMULATIONS AND DATA ANALYSIS 
2.1 Simulated galaxy images 

Charge Transfer Inefficiency has the greatest impact on 
small, faint objects that are far from the readout register (i.e. 


that have undergone a great number of transfers). To quan¬ 
tify the worst case scenario, we therefore simulate the small¬ 
est, faintest galaxy whose properties are likely to be mea¬ 
sured - with an exponential flux profile /(r) oc e”"^ whose 
broad wings (compared to a Gaussian or de Vaucouleurs 
profile) also make it more susceptible to CTI. To beat down 
shot noise, we simulate 10 ^ noisy image realisations for each 
measurement. We locate these galaxies 2048±0.5 pixels from 
both the serial readout register and the amplifier, uniformly 
randomising the sub-pixel centre to average out effects that 
depend on proximity to a pixel boundary. All our simulated 
galaxies have the same circularly symmetric proHle, follow¬ 
ing the observation by Rhodes et al. (20101 that this pro¬ 
duces the same mean result as randomly oriented elliptical 
galaxies with no preferred direction. 

We create the simulated images spatially oversampled 
by a factor 20 , convolve them with a similarly oversampled 
Point Spread Function (PSF), then resample them to the 
Hnal pixel scale. We use a preliminary PSF model and the 
0”l pixels of the Euclid VIS instrument, but our setup can 
easily be adapted to other instruments, e.g. ACS. To the 
image signal of ~1300 electrons, we add a uniform sky back¬ 
ground of 105 electrons, as expected for a 560 s VIS expo¬ 
sure, and Poisson photon noise to both the source and the 
background. After clocking and charge trailing (if it is being 
done; see Sect. 4.11, we then add additional readout noise, 
which follows a Gaussian distribution with a root mean 
square (rms) of 4.5 electrons, the nominal Euclid VIS value. 

In the absence of charge trailing, the final galaxies 
have mean S'/A’=11.35, and Full Width at Half Maximum 


(FWHM) size of 0'H8, as measured by SExtractor (Bertin 
|fc Arnouts||1996 1 . This size, the same as the PSF, at the 
small end of the range expected from Fig. 4 of [Massey et al.| 
(20131 makes our galaxies the most challenging in terms of 
CTI correction. Examples of input, degraded, and corrected 
images are shown in Fig. 

Separately, we perform a second suite of simulations, 
containing 10® realisations of a Euclid VIS PSF at S/N ~ 
200. The PSF simulations follow the above recipe, but skip 
the convolution of the PSF with an exponential disk. 


2.2 Image analysis 


On each of the sets of images (input, degraded, and cor¬ 
rected), we detect the sources using SExtractor. Moments 
of the brightness distribution and fluxes of the detected 
objects are measured using an IDL implementation of the 
RRG ( Rhodes, Refregier &: Groth]2001 1 shape measurement 
method. RRG is more robust than SExtractor for faint im¬ 
ages, combining Gaussian-weighted moments of the image 
I{0) to measure integrated source flux 


F=J w{e)i{e)d^0. 


( 1 ) 


where W{0) is a Gaussian weight function with standard 
deviation w, and the integral extends over 2.5w; the position 


y=J 6 ) 2 lF( 6 >)J( 6 »)d" 6 »; 

( 2 ) 

the size 


= Qii ± Q 22 ; 

( 3 ) 
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Figure 1. Example of three independent noise realisations of our simulated image of a faint galaxy; we used one million of these for each 
measurement. To the input images [upper left panel) the clocking is applied, degrading the image by both CTI and readout noise. The 
lower left panel shows the relative differences between the degraded and input images for a noise-free realisation [upper right panel). 
We can correct for the CTI, but not the readout noise (if present), by running the correction software with the trap model as used for 
clocking. If the correction model deviates from the clocking model, simulating our imperfect knowledge of the trap model, the correction 
will be imperfect. The lower right panel shows the relative difference image for a noise-free example, with the same scale as for the 
degraded image. We explore this imperfection; unnoticeable to the naked eye in noisy images, but crucial for Euclid'’^ performance. 


and the ellipticity 


HIGH PRECISION ELLIPTICITY 
MEASUREMENTS 


{61,62} = 


Qii — Q22 2Qi2 


Qii + Q22 ’ Qii + Q22 


( 4 ) 


where the second-order brightness moments are 


Qci 3=J ecepW{e)i{0)d^0, {a,p}e{i,2}. 


( 5 ) 


For measurements on stars, we chose a window size w = 0V75, 
the Euclid prescription for stars. For galaxies, we seek to re¬ 
produce the window functions used in weak lensing surveys. 


We adopt the radius of the SExtractor object (e.g. Lean 
thaud et al.|2007 1 that with w = 0”34 truncates more of the 


noise and thus returns more robust measurements. 

Note that we are measuring a raw galaxy ellipticity, a 
proxy for the (reduced) shear, in which we are actually in¬ 


terested (cf. Kitching et al. 2012 for a recent overview of 


the effects a cosmic shear measurement pipeline needs to 
address). A full shear measurement pipeline must also cor¬ 
rect ellipticity for convolution by the telescope’s PSF and 


calibrate it via a shear ‘responsivity factor’ (Kaiser, Squires 
|fc Broadhurst|[l995[ ). The first operation typically enlarges 
6 by a factor of ~ 1.6 and the second lowers it by about the 
same amount. Since this is already within the precision of 
other concerns, we shall ignore both conversions The abso¬ 
lute calibration of shear measurement with RRG may not be 
sufficiently accurate to be used on future surveys. However, 
it certainly has sufficient relative accuracy to measure small 
deviations in galaxy ellipticity when an image is perturbed. 


3.1 Measurements of a non-linear quantity 

A fundamental difficulty arises in our attempt to measure 
galaxy shapes to a very high precision, by averaging over a 
large number of images. Mathematically, the problem is that 
calculating ellipticity 6i directly from the moments and then 
taking the expectation value £{■) of all objects, i.e.: 


Cl = £ 


Qii — Q22 
Qii + Q22 


62 = £■ 


2Qi 


Qii + Q22 


( 6 ) 


means dividing one noisy quantity by another noisy quan¬ 
tity. Furthermore, the numerator and denominator are 
highly correlated. If the noise in each follows a Gaussian dis¬ 
tribution, and their expectation values are zero, the prob¬ 
ability density function of the ratio is a Lorentzian (also 
known as Cauchy) distribution. If the expectation values of 
the Gaussians are nonzero, as we expect, the ratio distribu¬ 
tion becomes a generalised Lorentzian, called the Marsaglia- 
Tin distribution ( Marsaglia|1965 2006 Tin[|1965 I. In either 
case, the ratio distribution has infinite second and first mo¬ 
ments, i.e. its variance - and even its expectation value - are 
undehned. Implications of this for shear measurement are 


et al. (20121; Kacprzak et al. (20121; Miller et al. (20131; 


discussed in detail by Melchior & Viola (20121; Refregier 


Viola, Kitching & Joachimi (20141. 


Therefore, we cannot simply average over ellipticity 
measurements for 10^ simulated images. The mean estima¬ 
tor (Eq. would not converge, but follow a random walk in 
which entries from the broad wings of the distribution pull 
the average up or down by an arbitrarily large amount. 
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3.2 “Delta method” (Taylor expansion) 
estimators for ellipticity 

As an alternative estimator, we employ what is called in 
statistics the ‘delta method’: a Taylor expansion of Eq. 
aronnd the expectation value of the denominator (e.g. 


Casella & Berger 2002 i. The expectation value of the ra¬ 
tio of two random variables X, T is thus approximated by: 


£{X/Yy, 


. £{X) 

‘ £{Y) 


C(X,Y) £{X)a 

£^iY) 

C{X,Yy 


'(Y) 


-b 


SHY) 


SHY) 

£{X)£[Y - £{Y)f 
£^{Y) 


(7) 


where £{X), u{X), ct^{X) denote the expectation value, 
standard deviation, and variance of X, and C(X,Y) its co- 
variance with Y. The zero-order term in Eq. Q is the often- 
used approximation £{X/Y) ~ £{X)/£{Y) that switches 
the ratio of the averages for the average of the ratio. We 
note that beginning from the first order there are two terms 
per order with opposite signs. Inserting Eq. © into Eq. 
the first-order estimator for the ellipticity reads in terms of 
the brightness distribution moments Qa/3 as follows: 


ei = 


S{Qi 1 —Q 22 ) 

S{Qii+Q22) 

O'^(Qll)—(T^(Q 22 ) SiQll—Q 22 )cr^{Qll+Q 22 ) 


62 = 


S‘^{Qi1+Q22 

S{ 2 Qi 2 ) 


SHQ11+Q22) 


( 8 ) 


S{Qii+Q 22 ) 

C{Qll, Qi2)+C{Qi2, Q22) 

f ^(Qii+Q22) 


+ 


S (2Qi2)a'^(Qii -bQ22) 


SHQi^+Q22) 


(9) 


with the corresponding uncertainties, likewise derived 
using the delta method (e.g. Casella fc Berger|2002 1: 


_2/. ^ cr^(Qii—Q22) 


He,) = 


S^{Q„+Q22) 

f(Qll —Q22) [o'^(Qll) —(J^(Q22)] 


-b - 


£’ 3 (gil-bQ 22 ) 

f ^(Qll — Q22)o'^(gil-bQ22) 


( 10 ) 


^’^(Qii -bg22) 


_2/. \ Cr^(Qii-bQ22) 


He2) = 


S^{Q„+Q22) 

f (Qii-bQ22) [C(Qii,Qi2)-bC(Qi2,Q22)] 

S^{Q„+Q22) ' 

^^^(gil+Q22)u^(gil-bQ22) 

shoHIy^) ■ 


( 11 ) 


3.3 Application to our simulations 

For our input galaxies, the combined effect of the first-order 
terms in eq. ([^ is ~ 10 %. Second-order contributions to the 
estimator are small, so we truncate after the first order. How¬ 
ever, because of the divergent moments of the Marsaglia-Tin 
distribution, the third and higher-order contributions to the 
Taylor series increase again. 

Nevertheless, while this delta-method estimator neither 


mitigates noise bias nor overcomes the infinite moments of 
the Marsaglia-Tin distribution at a fundamental level, it suf¬ 
ficiently suppresses the random walk behaviour for the pur¬ 
poses of this study, the averaging over noise realisations of 
the same object. We advocate re-casting the Euclid require¬ 
ments in terms of the Stokes parameters (Qii ± Q22,2Qi2', 
Viola, Kitching fc Joachimi]2014 1. These are the numerators 


and denominator of eq. ([6|) and are well-behaved Gaussians 
with finite first and second moments. 

The formal uncertainties on ellipticity we quote in the 
rest of this article are the standard errors (j{e,)/'/N given 
by eq. (10 b Our experimental setup of re-using the same 
simulated sources (computationally expensive due to the 
large numbers needed), our measurements will be intrinsi¬ 
cally correlated (Sect. |4.2| |. Hence the error bars we show 
overestimate the true uncertainties. 


4 THE EFFECTS OF FAST AND SLOW TRAPS 
4.1 How CTI is simulated 


The input images are degraded using a C implementation of 


the Massey et al. (20141 CTI model. During each pixel-to- 


pixel transfer, in a cloud of rig electrons, the number cap¬ 
tured is 


g(ne) = (1 - exp (^-anl j '^pi ’ 


( 12 ) 


where the sum is over different charge trap species with den¬ 
sity pi per pixel, and w is the full-well capacity. Parameter a 
controls the speed at which electrons are captured by traps 
within the physical volume of the charge cloud, which grows 
in a way determined by parameter /? . 

Release of electrons from charge traps is modelled by 
a simple exponential decay, with a fraction 1 — es¬ 

caping during each subsequent transfer. The characteristic 
release timescale n depends on the physical nature of the 
trap species and the operating temperature of the CCD. 

In this paper, we make the simplifying ‘volume-driven’ 
assumption that charge capture is instantaneous, so a = 0. 
Based on laboratory studies of an irradiated VIS CCD (de¬ 
tailed in Sect. 0 , we adopt a /3 = 0.58 baseline well fill, and 
end-of-life total density of one trap per pixel, 1. In our 
first, general tests, we investigate a single trap species and 
explore the consequences of different values of t. 


4.2 Iterative CTI correction 

The [Massey et J(] ( |2014[ ) code can also be used to ‘untrail’ 
the CTI. If required, we use niter = 5 iterations to attempt 
to correct the image (possibly with slightly different model 
parameters). Note that we perform this correction only after 
adding readout noise in the simulated images. 

Our main interest in this study is the impact of un¬ 
certainties in the trap model on the recovered estimate of 
an observable p (e.g. ellipticity). Therefore, we present our 
results in terms of differences between the estimators mea¬ 
sured for the corrected images, and the input values: 

Xpi = Pi 

,corrected input* (13) 

Because for each object of index i the noise in the measure¬ 
ments of Pi,corrected and pi,input are strongly correlated, they 
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Figure 2. The effect on measurements of galaxy flux F, astrometry y and morphology (size ellipticity ei) of charge traps of different 
species, i.e. release timescale r. The pixel-to-pixel transfer is assumed to be instantaneous. Filled (open) circles denote data for faint 
galaxies before (after) CTI mitigation. Crosses (plus signs) denote data for bright stars before (after) CTI mitigation. Data for CTI- 
corrected images are shown multiplied by a factor of 10 for Ay and Aei. Lines give the best fit following Eq. | |14[ |, with the coefficients 
listed in Tabled 
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Table 1. Parameters of fitting functions to illustrate the effect on measurements of galaxy fluxes F and Pg, astrometry y and mor¬ 
phology B?, e\ of charge traps of different species. In all cases, the measurements assume a density of one trap per pixel, and the 
astrophysical measurement is fitted as a function of the charge trap’s characteristic release time t as A -|- Da atan((logr — Dp)/Dai) + 
Gaexp((logT — Gp)^/2G^). Values after correction highlight the efficacy of CTI mitigation. 



A 

Da 

Dp 

D„ 

Ga 

Gp 

Gw 

Galaxy simulation: in degraded images, including readout noise [10 

X parameter] 




AF/ Ftrne 

-0.5367 ± 0.0098 

-0.3144 ± 0.0085 

6.199 ±0.044 

4.639 ± 0.260 

0.2116 ±0.0194 

49.53 ± 1.64 

41.54 ± 2.39 

Ay 

1.1098 ± 0.0014 

-0.5291 ± 0.0028 

8.392 ±0.080 

2.110 ± 0.234 

0.3061 ± 0.0185 

6.935 ±0.402 

7.083 ±0.210 


0.4226 ±0.0025 

-0.3857 ± 0.0038 

15.72 ± 0.18 

2.576 ±0.375 

1.0866 ±0.0448 

4.382 ±0.047 

3.779 ±0.160 

Aei 

0.5333 ±0.0016 

-0.3357 ± 0.0026 

16.28 ±0.22 

2.951 ± 0.326 

0.9901 ± 0.0203 

4.553 ±0.054 

4.132 ±0.081 

Galaxy simulation: after correction 

in software post-processing (perfect knowledge of charge traps) [100 x parameter] 


AF / Ftrue 

-0.5549 ± 0.0029 

0.0446 ± 0.0028 

129.6 ± 13.7 

26.00 ± 13.36 

0.1301 ± 0.0121 

73.47 ± 6.78 

56.84 ± 5.21 

Ay 

0.09582 ±0.01011 

0.0517 ±0.0111 

5.622 ±8.911 

2.227 ± 4.557 

0.0810 ±0.1170 

2.757 ± 5.369 

3.154 ± 2.784 

ARVRlae 

-2.3181 ± 0.0173 

0.4431 ± 0.0202 

75.90 ±25.02 

28.47 ± 11.03 

0.5471 ± 0.2294 

41.31 ± 16.09 

35.33 ±9.12 

Aei 

0.01383 ±0.0115 

0.0039 ± 0.0066 

12.30 ±20.49 

1.000 ± 0.000 

0.0982 ± 0.0274 

5.738 ±2.085 

5.353 ± 2.078 

Star simulation: in degraded images, including readout 

noise [100 X 

parameter] 




AF/ FtTue 

-2.2472 ±0.0239 

-1.4558 ± 0.0189 

107.5 ±0.3 

55.11 ± 0.95 

1.151 ± 0.047 

496.6 ±3.2 

343.6 ± 4.4 

Ay 

4.3532 ±0.0014 

-1.8608 ±0.0027 

173.1 ± 0.4 

29.20 ± 0.67 

5.0987 ± 0.0173 

67.20 ±0.20 

43.91 ± 0.22 

^RVRlas 

0.9489 ± 0.00098 

-6.434 ± 0.0095 

288.8 ±4.7 

18.71 ± 4.49 

20.237 ± 0.716 

94.42 ± 0.15 

50.20 ± 0.25 

Aei 

1.2336 ±0.0077 

-0.7941 ± 0.0086 

266.7 ± 2.4 

17.54 ± 3.90 

16.513 ±0.046 

94.87 ± 0.19 

52.57 ± 0.21 

Star simulation: after correction in 

software post-processing (perfect knowledge of charge traps) [lOOx 

parameter] 


AF/ Ftrue 

-0.0035 ±0.0002 

0.0027 ± 0.0003 

110.2 ± 10.5 

42.21 ± 20.02 

0.0006 ± 0.0271 

182.6 ±71.3 

3.5 ± 100.0 

Ay 

0.1504 ±0.00066 

0.0970 ± 0.0067 

12.46 ± 1.86 

2.731 ± 1.552 

0.0218 ±0.0034 

7.377 ± 1.024 

5.063 ±0.717 

ARVRla. 

-0.0163 ±0.0038 

-0.0182 ±0.0036 

1269 ± 33 

24.57 ± 47.63 

0.0198 ±0.0146 

50.83 ± 34.56 

37.95 ± 38.64 

Aei 

0.0012 ±0.0024 

0.0003 ± 0.0014 

2.26 ±50.92 

1.000 ± 0.000 

0.02668 ± 0.0061 

8.465 ± 1.800 

5.379 ± 1.647 


partially cancel out. Thus the actual uncertainty of each 
is lower than quoted. Moreover, because we re-use the same 
noise realisation in all our measurements (cases of different 
pi and Ti), these measurements are correlated as well. 


4.3 CTI as a function of trap timescale 


The impact of charge trapping is dependent on the defect 
responsible. Figure [^demonstrates the effect of charge trap 
species with different release times r on various scientific 
observables. To compute each datum (filled symbols), we 
simulate 10^ galaxies, add shot noise, add CTI trailing in 
the y direction (i.e. vertical in Fig.[^, only then add readout 
noise. Separately, we simulate 10® stars. Using eqs. ®“( 111, 


we measure mean values of photometry (top panel), astrom¬ 
etry (second panel) and morphology (size in the third, and 
ellipticity in the bottom panel). Our results confirm what 
Rhodes et al. (20101 found in a different context. 


Three trap regimes are apparent, for all observables. 
Very fast traps (r < 0.3 transfers) do not displace electrons 
far from the object; thus their effect on photometry is min¬ 
imal (top plot in Fig. j^. We observe significant relative 
changes in position, size, and ellipticity, forming a plateau at 
low T, because even if captured electrons are released after 
the shortest amount of time, some of them will be counted 
one pixel off their origin. This is probably an artifact: We 
expect the effect of traps with r < 0.1 to be different in 
an model simulating the transfer between the constituent 
electrodes of the physical pixels, rather than entire pixels. 

Very slow traps (r > 30 transfers) result in electrons 
being carried away over a long distance such that they can no 
longer be assigned to their original source image. Hence, they 
cause a charge loss compared to the CTI-free case. However, 


because charge is removed from nearly everywhere in the 
image, their impact on astrometry and morphology is small. 


The most interesting behaviour is seen in the transi¬ 
tional region, for traps with a characteristic release time of a 
few transfer times. If electrons re-emerge several pixels from 
their origin, they are close enough to be still associated with 
their source image, but yield the strongest distortions in size 
and ellipticity measurements. This produces local maxima 
in the lower two panels of Fig.[^ If these measurements are 
scientifically important, performance can - to some degree - 
be optimised by adjusting a CCD’s clock speed or operating 
temperature to move release times outside the most critical 
range 1 <t< 10 ( Murray et al.|2012 l. 


In the star simulations (crosses in Fig. [^ for degraded 
images, plus signs for CTI-corrected images), the CTI effects 
are generally smaller than for the faint galaxies, because the 
stars we simulate are brighter and thus experience less trail¬ 
ing relative to their signal. Still, we measure about the same 
spurious ellipticity Aei and even a slightly higher relative 
size bias AR^/Rtrue for the stars. The explanation is that 
the quadratic terms in the second-order moments (eq. [^ 
allow for larger contributions from the outskirts of the ob¬ 
ject, given the right circumstances. In particular, the wider 
window size w explains the differences between the galaxy 
and PSF simulations. Notably, the peak in the Aei(r) and 
AR^/Rlaeir) curves shifts from ~ 3 px for the galaxies to 
~ 9 px for the stars. Because the wider window function gives 
more weight to pixels away from the centroid, the photom¬ 
etry becomes more sensitive to slower traps. 


For a limited number of trap conhgurations, we have 
also tried varying the trap density or the number of trans¬ 
fers (i.e. object position on the CCD). In both cases, the 
dependence is linear. Overall, for all tested observables, the 
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measurements in the degraded images (Fig.[^ solid symbols) 
are well-fit by the empirical fitting function 


Table 2. The baseline trap model Ai. The model includes a base¬ 
line well fill power of /5o = 0.58. 


r) = p(^A + Da atan((logT - Dp)/Dw)-|- 

Gaexp((logr-Gp)V2G2)), (14) 

which combines an arc-tangent drop (“D”) and a Gaussian 
peak (“G”). The best fit-amplitudes {A, Ds, and Ga), posi¬ 
tions on the r axis (Dp and Gp) and widths (D„ and Gw), 
are listed in Table The same functional form provides a 
good match to the residuals after GTI correction, /'^®°‘d(p, r) 
(open symbols in Fig.[^. These residuals are caused by read¬ 
out noise, which is not subject to GTI trailing, but undergoes 
CTI correction (see Sect. I5.3T2 1. 


4.4 Predictive model for imperfect correction 

We set out to construct a predictive model of Ap, the 
CTI effect in an observable relative to the underlying true 
value (eq. |13[ ). There are two terms, the CTI degradation 
(eq. |14[ ), and a second term for the effect of the ‘inverse’ 
CTI correction allowing for a slightly imperfect CTI model: 




‘‘(p-l-Ap, r-b Ar). 


(15) 


Since CTI trailing perturbs an image by only a small 
amount, the correction acts on an almost identical image. 
Assuming the coefficients of eq. (141 to be constant, we get: 


A ^ ^degr + Ap, r + Ar) + (P, r), (16) 


where r^^{p,T) is approximately constant, and depends on 


the readout noise (see Section 5.31. We could expand this 


equation as a Taylor series, but the derivatives of / do not 
provide much further insight. 


Because eq. (121 is non-linear in the number of sig¬ 


nal electrons, our observation (Sect. 4.31 that the effects of 
CTI behave linearly in p is not a trivial result. Assuming 


this linearly in p, we can expand eq. (161 and factor out p. 


The combined effect of several trap species i with release 
timescales n and densities pi can then be written as: 


Baseline model 

i = l 

i = 2 

i = 3 

Trap density pi [px“^] 

0.02 

0.03 

0.95 

Release timescale [px] 

0.8 

3.5 

20.0 


data, with a full Euclid CTI model featuring multiple trap 
species (see Sect. We call this the M-I-5M experiment. 

Akin to Prod’homme et al. (20121 for Gaia, this study 
is useful in the larger context of the flow down of require- 


its imaging capabilities (Massey et al. 

20131 and instrument 

implementation (Cropper et al. 2013 

20141. In particular. 


cess will be determined both by its absolute instrumental 
performance and our knowledge about it. We now present 
the next step in the flow down: to what accuracy do we need 
to constrain the parameters of the Massey et al. (20141 CTI 
model? Future work will then determine which calibration 
observations are required to achieve this accuracy. 

While the final Euclid requirements remain to be con¬ 
firmed, we adopt the current values as discussed by |Cropper| 
|et al.| ( |2013| ). Foremost, the “CTI contribution to the PSF 
ellipticity shall be <1.1 x 10“‘^ per ellipticity component”. 

The Euclid VIS PSF model will bear an uncertainty due 
to CTI, that translates into an additional error on measured 
galaxy properties. For the bright stars (which have much 
higher S/N) tracing the PSF, Cropper et al.| (2013l quote a 
required knowledge of to a precision |cr(J?^)| <4 x lO”'^. 
We test this requirement with our second suite of simula¬ 
tions, containing 10® realisations of a Euclid VIS PSF at 
S/Nx200 (cf. Sec.j^l. 

In reality, CTI affects the charge transport in both CCD 
directions, serial and parallel. For the sake of simplicity, we 
only consider serial CTI, and thus underestimate the total 
charge trailing. There is no explicit photometric error budget 
allocated to CTI, while “ground data processing shall correct 
for the detection chain response to better than 0.7 % error 
in photometry in the nominal VIS science images”. 


Af\pi + Api,Ti + An) = {ri)+ 

i 

y] [pifH) - {pi + Ap/)f{Ti + An)], (17) 

in which we dropped the superscript of j-qj. g^ke of 
legibility. We are going to test this model in the remainder 
of this study, where we consider a mixture of three trap 


species. We find eq. (171 to correctly describe measurements 


of spurious ellipticity Aei, as well as the relative bias in 
source size AR? / BZ/rue ^-nd flux AF/Ftrue- 


5 EUCLID AS A CONCRETE EXAMPLE 
5.1 Context for this study 

To test the general prediction eq. ( |17[ ), we now evaluate the 
effect of imperfect CTI correction in simulations of Euclid 


5.2 CTI model for the Euclid VISual instrument 


Based on a suite of laboratory test data, we define a baseline 
model M of the most important CTI parameters {pi, Ti, Po). 
We degrade our set of 10^ simulated galaxies using M. The 
M + 5M experiment then consists of correcting the trailing 
in the degraded images with slight alterations to M. We 
investigate > 100 correction models M-l-dM, resulting in an 
impressive 1.4 x 10® simulated galaxies used in this study. 

Exposure to the radiation environment in space was 
simulated in the laboratory by irradiating a prototype of 
the e2v CCD273 to be used for Euclid VIS with a 10 MeV 


equivalent fluence of 4.8 x 10® protons/cm (Prod’homme 


et al.|[20l4 Verhoeve et al.|[2014 |. Characterisation experl 
ments were performed in nominal VIS conditions of 153 K 
temperature and a 70 kHz readout frequency. We refer to 
Appendix|^for further details on the experiments and data 
analysis. 

We emphasize that our results for ei pertain to 
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corrected images, S/N = 11 galaxies 


corrected images, S/N=200 stars 
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Figure 3. Residual spurious ellipticity and source size induced by CTI after correction, as a function of the rms readout noise. Left 
plot: Spurious ellipticity Aei in the faintest galaxies to be analysed for Euclid, at the end of the mission, and furthest from the readout 
register. Right plot: Bias in source size in bright stars, at the end of the mission, and furthest from the readout register. 

Lines show the best quadratic (cubic) fits. Shaded regions indicate the Euclid VIS requirements. Vertical grey lines mark the nominal 
rms readout noise of 4.5 electrons. 


faint and small galaxies, with an exponential disk profile 
(vz. Sect. 2.11, and placed at the maximum distance from 
the readout register {y = 2051 transfers). Furthermore, we 
assume the level of radiation damage expected at the end 
of Euclid’s six year mission. Because CTI trailing increases 
roughly linearly with time in orbit (cf. [Massey et al.|20'l4 1 , 
the CTI experienced by the typical faintest galaxy (i.e. at 
half the maximum distance to the readout register and three 
years into the mission), will be smaller by a factor of 4 com¬ 
pared to the results quoted below. 

Where not stated otherwise the nominal Euclid VIS rms 
readout noise of 4.5 electrons was used. Table [2] summarises 
the baseline model M that was constructed based on these 
analyses. The default well fill power is /?o = 0.58. Slow traps 
with ra =20 clock cycles and ps = 0.95 dominate our baseline 
model, with small fractions of medium-fast (t 2 = 3.5, p 2 = 
0.03) and fast (n = 0.8, pi = 0.02) traps. Figure A3 shows 
how trails change with changing trap parameters. 


5.3 Readout noise impedes perfect CTI correction 

5.3.1 Not quite there yet: the zeropoint 

First, we consider the ellipticities measured in the degraded 
and corrected images, applying the same baseline model 
in the degradation and correction steps. The reasons why 
this experiment does not retrieve the same corrected ellip¬ 
ticity CcoTT as input ellipticity ein are the Poissonian im¬ 
age noise and Gaussian readout noise. We quantify this in 
terms of spurious ellipticity Ae = ecorr — Cin, and shall re¬ 
fer to it as the zeropoint of the M + 5M experiment. The 
spurious ellipticity in the serial direction is Ze^ = Aei = 
—0.00118 ± 0.00060. Thus, this experiment on worst-case 
galaxies using the current software exceeds the Euclid re¬ 
quirement of |Aea| < 1.1 X 10“'* by a factor of ~ 10. With 
respect to the degraded image 99.68 % of the CTI-induced 
ellipticity are being corrected. Virtually the same zeropoint, 
Aei =—0.00118±0.00058, is predicted by adding the contri¬ 
butions of the three species from single-species runs based 


on the full 10' galaxies. We point out that these results on 
the faintest galaxies furthest from the readout register have 


been obtained using non-flight readout electronics (cf. Short 
jet al.|2014| ). 

From our simulation of 10® bright (S/N « 200) stars, 
we measure the residual bias in source size after CTI 
correction of = AR^/R^rue = (—0.00112 ± 0.00030), in 
moderate excess of the requirement | AR^/Rj^uej <4 x 10“'^. 
While the S/N of the star simulations is selected to represent 
the typical Euclid VIS PSF tracers, the same arguments 
of longest distance from the readout register and non-flight 
status of the electronics apply. 


5.3.2 The effect of readout noise 

In Fig. we explore the effect of varying the rms readout 
noise in our simulations about the nominal value of 4.5 elec¬ 


trons (grey lines) discussed in Sect. 5.3.1 We continue to use 
the baseline trap model for both degradation and correction. 
For the rms readout noise, a range of values between 0 and 
15 electrons was assumed. For the faint galaxies (Fig.|^ left 
plot), we find Aei to increase with readout noise in a way 
well described by a second-order polynomial. A similar, cu¬ 
bic fit can be found for AR^/Rtrue measured from the star 
simulations (Fig.|^ right plot), but with a hint towards sat¬ 
uration in the highest tested readout noise level. 

The most important result from Fig. is that in ab¬ 
sence of readout noise, if the correction assumes the cor¬ 
rect trap model M, it removes the CTI perfectly, with 
Aei = (0.3±5.9) x IQ-^ and ARVRLe = (0.0±2.8) x 10“'‘. 
The quoted uncertainties are determined by the N = 10^ 
(10®) galaxy images we simulated. We conclude that the 
combination of our simulations and correction code pass this 
crucial sanity check. If the rms readout noise is < 3 electrons 
(< 0.5 electrons), the spurious ellipticity (the relative size 
bias) stays within Euclid requirements. 
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Figure 4. Sensitivity of the CTI-induced spurious ellipticity Aei 
(upper plot) and the relative spurious source size 
(lower plot) to the well fill power /J. At the default value of 
/3 = 0.58 (vertical grey line), the measurements deviate from zero 
due to readout noise, as indicated by arrows. The shaded region 
around the measurements indicate the Euclid requirement ranges 
as a visual aid. Solid and dashed lines display quadratic (linear) 
fits to the measured Aei(/3) and AiJ^(/3)/i?^^^g, respectively. We 
study the worst affected objects (at the end of the mission and 
furthest from the readout register) and the faintest Euclid galax¬ 
ies. This plot also assumes CTI is calibrated from charge injection 
lines at full well capacity only. This will not be the case. 


5.4 Sensitivity to imperfect CTI modelling 

5.4-1 Morphology biases as a function of well fill power, 
and determining tolerance ranges 

Now that we have assessed the performance of the correction 
using the same CTI model as for the degradation (given the 
specifications of our simulations), we turn to the M + 5M 
experiment for determining the sensitivities to imperfections 
in the CTI model. To this end, we assume the zeropoint 
offset Ze 


LJei (or Zfj 2 ) of Sect. 5.3.1 to be corrected, and ‘shift’ 
the requirement range to be centred on it (see, e.g.. Fig. [^. 

Figure shows the M -|- 5M experiment for the well 
fill power f). If the degraded images are corrected with the 
baseline /3o = 0.58, we retrieve the zeropoint measurement 
from Sect. For the M -I- dM experiment, we corrected 

the degraded images with slightly different well fill powers 


0.56^/35(0.60. The upper plot in Fig.fflshows the resulting 
Aei in galaxies, and the lower plot Ai?/i3true in stars. We 
find a strong dependence of both the spurious serial elliptic¬ 
ity Aei and AR^/R^,.„g on A/3 = /3—/?o. 

In order to determine a tolerance range with respect to 
a CTI model parameter ^ with baseline value (here, the 
well fill power (3), we fit the measured bias A 77 (e.g. Aei, cf. 
eq. 131 as a function of A^ = ^—^ 0 . By assuming a polynomial 




(18) 


1=1 


of low order J, we perform a Taylor expansion around ^o- 
In eq. 18 Z,, is the zeropoint (Sect. 5.3.11 to which we 


have shifted our requirement margin. The coefficients Oj 
are determined using the IDL singular value decomposition 
least-square fitting routine SVDFIT. For consistency, our fits 
include Zr, as the zeroth order. In Fig. the best-fitting 
quadratic (linear) fits to Aei /Rtrue) S'l'® shown as a 

solid and dashed line, respectively. 

In both plots, the data stick tightly to the best-fitting 
lines, given the measurement uncertainties. If the measure¬ 
ments were uncorrelated, this would be a suspiciously or¬ 
derly trend. However, as already pointed out in Sect. |3.3[ 
we re-use the same 10 ^ simulations with the same peaks and 
troughs in the noise in all data points shown in Figs. |4]to|^ 
Hence, we do not expect to see data points to deviate from 
the regression lines to the degree their quoted uncertainties 
would indicate. As a consequence, we do not make use of 
the Xred ^ 1 our fits commonly yield for any interpretation. 

Because the interpretation of the reduced is tainted 
by the correlation between our data points, we use an alter¬ 
native criterion to decide the degree J of the polynomial: If 
the uncertainty returned by SVDFIT allows for a coefficient 
Oj =0, we do not consider this or higher terms. For the pan¬ 
els of Fig.|^ this procedure yields J = 2 (J = 1). The different 
signs of the slopes are expected because Rd appears in the 
denominator of eq. Q. 

Given a requirement Arjreq, e.g. Aei,req = 1.1 x 10“’^, the 
parametric form (eq. |18| l of the sensitivity curves allows us 
to derive tolerance ranges to changes in the trap parameters. 
Assuming the zeropoint (the bias at the correct value of ^) 
to be accounted for, we find the limits of the tolerance range 
as the solutions A^toi of 




= ATjr, 


(19) 


with the smallest values of | A^| on either sides to A^ = 0 . 
Using, eq. (191, we obtain A/3toi = ±(6.31±0.07) x 10“® from 
the requirement on the spurious ellipticity Aei < 1.1 x 10“’^, 
for which the quadratic term is small. From the requirement 
on the relative size bias Ai?^/i?t,.ue < 4 x 10“'* we obtain 
A/3toi = ±(4.78±0.05) x 10“'*. In other words, the ellipticity 
sets the more stringent requirement, and we need to be able 
to constrain /3 to an accuracy of at least (6.31 ± 0.07) x 
10“® in absolute terms. This analysis assumes calibration 
by a single charge injection line at full well capacity, such 


that eq. (121 needs to be extrapolated to lower signal levels. 


We acknowledge that Euclid planning has already adopted 
using also faint charge injection lines, lessening the need to 
extrapolate. 
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Figure 5. Sensitivity of the CTI-induced spurious ellipticity Aei in faint galaxies (upper panel) and relative bias in source size AR? / 
(lower panel) in bright stars to a relative change in trap densities [pi + Api)/pi. Different symbols and line styles denote which to which 
of the trap species a change in density was applied: The slow traps: (pi, upward and dashed line); the medium traps: (p 2 , downward 
triangles and long-dashed line); both of them: (pi, p 2 , diamonds and triple dot-dashed line); the fast traps (p 3 , squares and dot-dashed 
line); all: (pi, p 2 , p 3 , circles and dotted line). The various broken lines show the best-fit representation of the measurements as given 
by the sensitivity model (Eq. |20[ |. Like in Fig. the grey shaded area indicates the Euclid VIS requirement range. We study the worst 
affected objects (at the end of the mission and furthest from the readout register) and the faintest Euclid galaxies. 
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5.4-S Ellipticity bias as a function of trap density 

We now analyse the sensitivity of Ae^ towards changes in 
the trap densities. Figure]^ shows the M + 5M experiment 
for one or more of the trap densities pi of the baseline model. 
The upper panel of Fig. presents the spurions ellipticity 
Aei for five different branches of the experiment. In each of 
the branches, we modify the densities pi of one or several 
of the trap species. For example, the upward triangles in 
Fig. [^denote that the correction model applied to the images 
degraded with the baseline model used a density of the fast 
trap species pi+Api, tested at several values of Api with 
0.9 ^ 1 + Api/pi 1.1. The densities of the other species 
are kept to their baseline values in this case. The other four 
branches modify p2 (downward triangles); pa (squares); pi 
and p 2 (diamonds); and all three trap species (circles). 

Because a value of Api = 0 reproduces the baseline 
model in all branches, all of them recover the zeropoint 


measurement of Aei there (cf. Sect. 5.3.11. Noticing that 


Sdegr,! — ein,i <0 for the degraded images relative to the in¬ 
put images, we explain the more negative Aei for Api < 0 
as the effect of undercorrecting the CTI. This applies to 
all branches of the experiment. Likewise, with increasing 
Api > 0 , the residual undercorrection at the zeropoint de¬ 
creases. Eventually, with even higher k> 1, we overcorrect 
the CTI and measure Aei >0. 

Over the range of 0.9^ H-Api/pi ^ 1.1 we tested, Aei 
responds linearly to a change in the densities. Indeed, our 


model (eq. 171, which is linear in the pi and additive in the 
effects of the different trap species, provides an excellent de¬ 
scription of the measured data, both for Aei and AR^ / 
(Fig.d lower panel). The lines in Fig. [^denote the model 
prediction from a simplified version of eq. flTf, 

A/P'(p, + Ap,) = ^ P.r=‘‘'(r,) + ^(P. + Ap,)/(ri). (20) 


In eq. (20 1 , we assumed the ri be correct, i.e. Ari = 0. 


Next, we compute the tolerance Api,toi/p by which, 
for each branch of the experiment, we might deviate from 
the correct trap model and still recover the zeropoint 
within the Euclid requirements of |Aec,,req| < 1.1 x 10”“^, 
resp. I Ai?req/i?true I < 4x 10 “"^. Again, we calculate these tol¬ 


erances about the zeropoints Z — pif^^’^"^(ri) (cf. eq. 20 1 , 

that we found to exceed the requirements in Sect. |5.3Tlj but 
assume to be corrected for in this experiment. 

In accordance with the linearity in Api, applying the 
Taylor expansion recipe of Sect. |5.4.H we find the data 
in F ig. [5| to be well represented by first-order polynomials 
(eq. |18^T The results for Api^toi/p we obtain from eq. ( [T^ 
are summarised in Table For all species, the constraints 
from Aei for faint galaxies are tighter than the ones from 
AR? /AR^^^g for bright stars. 

Only considering the fast traps, pi can change by 0.84± 
0.33% and still be within Euclid VIS requirements, given the 
measured zeropoint has been corrected for. While a tolerance 
of 0.39 ± 0.06% is found for p 2 , the slow traps put a much 
tighter tolerance of 0.0303±0.0007% on the density pa. This 
is expected because slow traps amount to 95% of all baseline 
model traps (Table[^. Varying the density of all trap species 
in unison, we measure a tolerance of 0.0272 ± 0.0005%. 

Computing the weighted mean of the At = 0 intercepts 
in Fig. 1^ we derive better constraints on the zeropoints: 


= Aei = —0.00117 ± 0.00008 for the faint galaxies, and 
=Ai?^/i?true = —0.00112 ± 0.00004 for the bright stars. 


5 . 4.3 Ellipticity bias as a function of trap release time 

Figure shows the M ± 5M experiment for one or more 
of the release timescales ri of the trap model. The upper 
panel of Fig. [^presents the spurious ellipticity Aei for five 
different branches of the experiment. In each of the branches, 
we modify the release timescales Ti of one or several of the 
trap species by multiplying it with a factor {n + Ari)/ri. 

As in Fig.[^ the upward triangles in Fig. [^denote that 
the correction model applied to the images degraded with 
the baseline model used a density of ri ± Ari for the fast 
trap species. The release timescales of the other species are 
kept to their baseline values in this case. The other four 
branches modify T 2 (downward triangles); T 3 (squares); ri 
and T 2 (diamonds); and all three trap species (circles). 

Because a value of At = 0 reproduces the baseline 
model in all branches, all of them recover the zeropoint 
measurement of Aei there. The three trap species dif¬ 
fer in how the Aei they induce varies as a function of 
Ari. One the one hand, for ri, we observe more negative 
Aei for (t; ± Ari)/ri < 1 , and less negative values for 
(n ± ATi)/Ti > 1, with a null at (Ti ± Ari)/ri « 1.5. On 
the other hand, with the slow traps (ra), we find Aei > 0 
for (Ti + ATi)/ri <0.99, and more negative values than the 
zeropoint for (Ti + Ar^)/ri > 1. The curve of Aei('^''" 2 ) shows 
a maximum at (ri ± Ari)/Ti« 0 . 8 , with a weak dependence 
on 0.7<(ri ± Ari)/ri<l.l. 

Key to understanding the spurious ellipticity as a func¬ 
tion of the Ti is the dependence of Aei(T) for a single trap 
species that we presented in Fig.[^ and expressed by the em¬ 
pirical fitting function /e^ (r) (Eq. |14[ ) with the parameters 
quoted in Table While the correction algorithm effectively 
removes the trailing when the true Ti is used, the residual 
of the correction will depend on the difference between the 
Aca for Ti and for the timescale (n ± ATi)/Ti actually used 
in the correction. This dependence is captured by the pre¬ 


dictive model (Eq. 171, which simplifies for the situation in 
Fig. [^(Ap, =0) to 


A/’^''(r* ± Ari) = Z ± [f(Ti) - f(Ti + An)] , 


( 21 ) 


with Z = (lines in Fig. |^. In the branches 

modifying ti and/or T 2 , but not ra, the measurements over 
the whole range of 0.5 ^ (Ti ± ATi)fTi ^ 1.6 agree with the 
empirical model within their uncertainties. If ra is varied, 
Eq. ( |21[ ) overestimates jAeij significantly for j An] >O.OSn. 
We discuss a possible explanation in Sect. Our empirical 
model provides a natural explanation for the maximum in 
Aei(T2): Because T2 = 3.5 is located near the peak in /ei(T), 
assuming (n ± An)/n ^ 0.8 for correction means using a 
release time regime where Aei(r) is still rising instead of 
falling. The correction software accounts for this; hence the 
spurious ellipticity from using the wrong release time scale 
shows the same maximum as /ei(T). 

Because T 2 is not located very closely to the peak in 
AR^ARlUr) (cf. Fig. [^, we do not see an extremum in 
the lower panel of Fig. which shows the sensitivity of the 
size bias in bright stars to variations in the Ti. 

In order to compute the tolerances Artoi towards 
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Figure 6. Sensitivity of the CTI-induced spurious ellipticity Aei in faint galaxies (upper panel) and relative bias in source size AR? / 
(lower panel) in bright stars to a relative change in release timescales (r^ + ATi)/Ti. Different symbols and line styles denote to which of 
the trap species a change in release timescale was applied: The slow traps: (ri, upward triangles and dashed line); the medium traps: (r 2 , 
downward triangles and long-dashed line); both of them: (ri, T 2 , diamonds and triple dot-dashed line); the fast traps (t 3 , squares and 
dot-dashed line); all: (ri, r 2 , ra, circles and dotted line). The various broken lines show the best-fit representation of the measurements 
as given by the empiric sensitivity model (eq. |21| l. Like in Fig.j^ the grey shaded area indicates the VIS Euclid requirement range. We 
study the worst affected objects (at the end of the mission and furthest from the readout register) and the faintest Euclid galaxies. 


© 2014 RAS, MNRAS 000, 
















How well can CTI be corrected? 13 


Table 3. Tolerances for changes in the trap model parameters ^ (column 1), derived from polynomial fits to the sensitivity curves, and 
taking into account the Euclid VIS requirements (eqs. |18[|19[ |. Here, we assume no residual biases when using the correct trap model. For 
three observables 77 , we show pairs of tolerances, for < 0 and > 0 each. Columns 2 and 3 report tolerances based on the ellipticity 
bias Aei, columns 4 and 5 those from the relative size bias AR ^/and columns 6 and 7 those from the relative flux bias AF/Ftrue. 


branch ^ 

10 " ACS” 

lO^AC^S 

10 "AC”i” 

lO^AC^S 

10" ACS,” 

10'‘AC“S 


77 = Aei 

77 = Aei 

7 = AR 2 /ij 2 ^^^ 

77 = AR 2 /ij 2 ^^^ 

77 = AF/Ftrue 

77 = AF/Ftrue 


galaxies 

galaxies 

stars 

stars 

galaxies 

galaxies 

0 

-0.631 ± 0.007 

0.631 ± 0.007 

-4.78 ±0.05 

4.78 ± 0.05 

-61.5 ± 0.3 

60.5 ± 0.3 

pi 



-12501«°o 

125011°°° 

— 

— 

P2 

-39l^ 

391° 

-19ll(° 

19llii^ 

— 

— 

P3 

Q r\q + 0.06 

3-03l°;°^ 

-5.91 ± 0.03 

5.91 ± 0.03 

-267.5 ± 1.6 

267.5 ± 1.6 

Pi,2 

-2612 

26l° 

-166112 

166111 

— 

— 

Pi,2,3 

-2.72 ±0.05 

2.72 ± 0.05 

-5.71 ± 0.03 

5.71 ± 0.03 

-262.8 ± 1.6 

262.8 ± 1.6 

TI 

-193119 

19312° 

-13101120 

1310111° 

< -10000 

> 10000 

T2 

-300l9°o 

270l7o° 

-270l” 

270l°° 

< -10000 

> 10000 

T3 

-4.00 ±0.04 

4.00 ± 0.04 

-11.30 ±0.05 

11.31 ± 0.05 

-1574+11 

2320+11° 

Ti,2 



-220l°° 

230l°° 

< -10000 

> 10000 

Ti,2,3 

-4.03 ±0.04 

4.04 ± 0.04 

-11.69 ±0.05 

11.68 ±0.05 

-145411° 

20201™ 

ti,2,3, pi, 2 , 3 , first pixel matched 

-16.071°;” 

16.091°;” 

-16.17 ± 0.09 

16.21 ± 0.09 

-262.5 ±0.7 

263.0 ± 0.7 


changes in the release timescales, we again employ a polyno¬ 
mial fit (eq. 191. Evidently, the tolerances differ substantially 
between the r^, again with the narrower tolerance intervals 
from Aei than from AR ^/Only for At 2 with its 
extreme point for Aei near the baseline value, we find simi¬ 
lar tolerances In both cases. However, even for the rare trap 
species ti, the tolerance Is only Ari,toi = (1.93 ± 0.23) %. 
One needs to know the release timescale of the slow trap 
species to an accuracy of (0.0400 ± 0.0004) % to be able to 
correct it within Euclid VIS requirements. We find the same 
tolerance if all timescales are varied in unison. 


5.4-4 Combinations of timescales and densities yielding 
the same first trail pixel flux 

Considering how trap parameters are constrained practically 
from Extended Pixel Edge Response (EPER) and First Pixel 
Response (FPR) data, it is instructive to consider combina¬ 
tions of trap release timescales Ti and densities pi that yield 
the same number of electrons in the first pixel of the trail as 
the baseline model. This is interesting because given realis¬ 
tic conditions, the first pixel of the trail will have the largest 
signal-to-noise ratio and will be most easily constrained. We 
thus perform an initial exploration of the parameter degen¬ 
eracies. In our “first pixel matched” models, the effect of a 
given change in r on the first trail pixel needs to be compen¬ 
sated by a change in p. Because a larger (smaller) r means 
more (less) charge close to the original pixel, the compensa¬ 
tion requires Ap^ < 0 for Art < 1 and Ap; > 1 for Axi > 1. 
Only in the branches where we vary T 3 or all timescales to¬ 
gether, we And the Api to differ noticeably from unity. For 
the latter two, they populate a range between Api = 0.745 
for Ari = 0.7 to Api = 1.333 for Ar = 1.4. 

Figure shows the M -I- dM experiment for all Ti and 
pi (large symbols). Small symbols depict the alteration to 
Ti, but with the pi kept fixed, i.e. the same measurement as 
the open circles in Fig. Compared to these, Aei in faint 
galaxies (upper panel) is of opposite sign in the “first pixel 
matched” case, relative to the zeropoint. This can be under¬ 
stood as an effect of our baseline trap mix being dominated 


by slow traps, for which a small increase in t leads to less 
CTI-lnduced ellipticity. The simultaneous increase in trap 
density effects more CTI-induced ellipticity, and this is the 
larger of the two terms, such that a change in sign ensues. 
The same holds for AR^/Rt^ue in bright stars (lower panel 
of Fig. [^, but with inverted slopes compared to Aei. 

Again using eq. (191, we compute the tolerance range 
for the changes to the Ti in the “first pixel matched” case. 
(The respective changes to the pi are determined by the first 
pixel constraint.) Modifying all release time scales, we arrive 
at Artoi = 0.16 %. (Table [^. This tolerance is wider than 
the 0.04% for Aei when only the Ti are varied, again due to 
the different signs arising from variations to T 3 and ps- By 
coincidence, we also arrive at Artoi =0.16 % when repeating 
that test with the size bias measured in bright stars. 

The black solid line in Fig.[^shows the predictive model 
(eq. 171, taking into account the combined effect of the Art 
and Api, giving the same first pixel flux. Both in the ri-only 
(dotted line) and “first pixel matched” cases it matches the 
measurements only within a few percent from A = 1. Cru¬ 
cially, this mismatch only occurs for Aei in faint galaxies, 
but not for AR^/Rt^og in bright stars. 

We explain this discrepancy with the uncertainties with 
which our measurements and modelling (Fig.|^ describe the 
underlying function /ei(r). The range 20<r<100 is where 
the fitting function Eq. (141 deviates most from the obser¬ 
vations in Fig. The CTI correction effectively removes 
almost all CTI effects on photometry and morphology, leav¬ 
ing the residuals presented in Figs. [5]to[^ at least one order 
of magnitude smaller than the scales of the uncorrected CTI 
effects. Hence, a relatively small uncertainty in fir) causes 
a large mismatch with the data. 

The cause of the uncertainty in the parameters of 
Eq. (14|, shown in Table is twofold: First, there is uncer¬ 
tainty in the fit as such. Second, there is uncertainty due to 
the finite sampling of the AeaH) AR.ei('r) curves. Run¬ 
ning a denser grid in r can remove the latter, but the former 
might be ultimately limited by our choice of the function 
(Eq. |14[ ), which is empirically motivated, not physically. We 
further discuss the limits of the predictive model in Sect. 
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Figure 7. The same as Fig. but for Ar^ < 0 combinations of 
timescales and densities pi that yield the same count rate in the 
first trail pixel as the baseline model. All trap species are modified 
in unison (large symbols and solid line). For comparison, small 
symbols and the dotted line repeat the result from Fig. where 
only the r; were modified, not the pi. (Notice the different scale 
of the ordinates.) The lines show the predictive models (Eg. |17[ |. 
We study the worst affected objects (end of the mission, furthest 
from the readout register) and the faintest Euclid galaxies. 


Figure 8. Relative bias in RRG flux with respect to the true 
input flux, as a function of readout noise (upper panel) and well 
fill power 0 (lower panel). Solid lines give the best-fit polynomial 
models. The grey-shaded Euclid requirement range is centred on 
zero for the readout noise plot, and on the zeropoint correspond¬ 
ing to the default readout noise for the 0 plot. Measurement un¬ 
certainties are shown, but very small. We study the worst affected 
objects (end of the mission, furthest from the readout register) 
and the faintest Euclid galaxies. 


5.5 Residual flux errors after imperfect CTI 
correction 

5.5. 1 Flux bias as a function of readout noise 

Given the default rms readout noise of 4.5 electrons, we 
measure a flux bias AF^si = AF/Ftiue relative to the 
true flux Ftrue in the input faint galaxy simulations of 
(—1.980 ± 0.012) % after CTI correction, corresponding to 
92.9 % of the CTI-induced flux bias being corrected. The up¬ 
per panel of Fig. shows the relative flux biases before and 
after correction as a function of rms readout noise. With¬ 
out readout noise, the flux bias can be corrected perfectly 
(Afj-ei = (0.002 ± 0.012) X 10“^ after correction). With in¬ 
creasing readout noise, the flux bias deteriorates, in a way 
that can be fitted with a cubic polynomial in terms of read¬ 
out noise. Comparing to the degraded images, we notice that 
the correction software applies same amount of correction. 


independent of the readout noise. Because the mitigation al¬ 
gorithm in its current form does not include a readout noise 
model, this confirms our expectations. 

We show the Euclid requirement on photometric accu¬ 
racy as the grey-shaded area in Fig.(upper panel), centred 
on zero. The nominal readout noise case exceeds the require¬ 
ment of <0.7 % photometric uncertainty for the faintest, 
worst-affected galaxies we study. However, the CTI-induced 
bias affects all VIS images, and would thus be calibrated out. 
The Euclid flux requirement can be understood as pertain¬ 
ing to uncertainties, not biases in the photometric calibra¬ 
tion. The uncertainty of the flux bias, 0.0012 % then makes 
only a tiny contribution to the photometric error budget. We 
now go on to study the sensitivity of the flux bias towards 
changes in the trap model. 
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5.5.2 Flux bias as a function of well fill power fi 

The lower panel of Fig. shows how a change in well fill 
power P alters the flux bias. If we correct the degraded im¬ 
ages nsing a p > fio, the model accounts for less CTI in 
small charge packages, i.e. less CTI in the image’s wings that 
are crncial for both photometry and morphology (cf. fig. 
Hence, a. P> Po leads to an undercorrection relative to the 
flux bias zeropoint Zf (Sect. I .5.1[ ), while for P—Po < —0.017, 
the zero line is crossed and overcorrection occurs. 

Although AFrei{P) in Fig. appears linear, using the 
criterion based on significant components (Sect. 5.4.21, a 


quadratic is preferred, indicated by the solid line. Us¬ 
ing eq. (191, we compute the tolerance range in P given 


AF’rei(/3toi) = 0.007, centred on Zf. Towards smaller well fill 
powers, we find A/3™J" = — (6.15±0.03) x 10~^, while towards 
larger p, we find A/3™^ = (6.05 ± 0.03) x 10“®. Compared 
to the constraints on the knowledge of P from Aei derived 
in Sect. |5.4.1[ these margins are ~ 100 times wider. 


5.5.3 Flux bias as a function of trap densities 

The upper plot of Fig.[^shows the flux bias AF^ei in depen¬ 
dence of a change Api to the densities pi in the correction 
model, in analogy to Sect. |5.4.2| Unless the density of the 
dominant trap species ps is modified, we measure only in- 
signihcant departures from the zeropoint Zf. Given the high 
accuracy of the flux measurements, these are still significant 
measurements, but they are negligible with respect to the 
Euclid requirement on flux. If all pi are varied in unison, 
the effect on ATrei is largest. A linear model using Eq. (181 
yields a tolerance of Apl°^jpi =±2.628±0.016 %, wider than 
the tolerances for Aei or AB? /(Table 1^. The lines in 
the upper plot of Fig.j^show that the model (eq. |20| | matches 
our measurements well over the range in Api we tested. 


5 . 5.4 Flux bias as a function of release timescales 

The lower plot of Fig. shows the flux bias AErei in de¬ 
pendence of a change An in the correction model, like in 


Sect. 5.4.3 As for varying the pi (Sect. 5.5.31, a change to 
only the fast and/or the medium traps yields only small 
departures from the zeropoint such that we can bundle to¬ 
gether all trap species for deriving a tolerance range. The re¬ 
spective measurements (black circles in Fig.[^ show a steep 
slope at An < 0 that flattens out to An > 0. This can be 
explained given the saturation of AErei found at large r in 
Fig. 2 and is confirmed by our model (eq. 21 dotted line in 
Fig. 91. Our prediction is offset from the measurement due 
to uncertainties in the modelling, but the slopes agree well. 


Although polynomial fits using eq. (181 warrant cubic 
terms in both cases, AF,.ei{n + An) is much straighter in 
the “first pixel matched” case where also the pi are altered 
(star symbols in Fig.[^ cf. Sect. 5.4.41. The reason is that the 
slopes of AF.rei{pi+Api) and AFj.ei{n+An) have the same 
sign and do not partially cancel each other out, as is the case 


for Aei(pi-I-Api) and Aei(ri-I-Ari). Again, eq. (171 succeeds 
in predicting the measurements, despite offsets that are sig- 
nihcant given the small uncertainties but small in terms of 
AFrei in the uncorrected images. 

Using the cubic fits, we find the following wide tolerance 


ranges (eq. 191 = 15.7 ± 0.2 % and = 



- 0.01 


- 0.02 


- 0.03 
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Figure 9. Upper panel: The same as Fig. but showing the 
sensitivity of the measured flux bias AF/Ftme as a function of 
the relative change in trap densities pi. Lower panel: The same 
as Fig.[^ but showing the flux bias AF/Ftme as a function of the 
relative change in trap densities r^. Star symbols and the solid 
line denote the “first pixel matched” model for all trap species. 
The lines in both panels show the predictive model (eq. [13 . We 
study the worst affected objects (end of the mission, furthest from 
the readout register) and the faintest Euclid galaxies. 


23.2/jJ g %. In the “first pixel matched”, case the intervals 
are considerably tighter, due to the contribution from the 
change in densities, with Ar/’^in/ri = 2.625 ± 0.007 % and 
^U^max/'Ti = 2.630±0.007 %. Again, the strictest constraints 
come from the ellipticity component Aei. 


6 DISCUSSION 

6.1 Limits of the predictive model 

We measured tolerance ranges for changes in the pi and n 
given the Euclid VIS requirements, and presented a model 
(Eq. |17[ ) capable of predicting these results based on the 
Ap(t) curves (e.g. Aei(r), Fig. [^, that are less expensive 
to obtain in terms of CPU time. However, as can be seen in 
particular in Fig.[^ there is a mismatch between predictions 
and measurements for rs, the most common baseline model 
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trap species. As discussed in Sect. |5.4.4[ this is caused by 
the finite sampling and the empirical nature of eq. (14l. 
Unfortunately, /(r) and will likely depend non- 


trivially on the source profile. Moreover, Eq. ( |17[ ), if applied 
to ellipticity, treats it as additive. Where this approximation 
breaks down, i.e. when values that are not ^ 1 are involved, 
the correct additional formula (e.g. |Schneidei]|2006[ ) must 
be used. This applies to CTI-induced ellipticity as well as to 
large intrinsic or shear components. 

We tested that the dependence on /3 (Fig. can be 
included in the model as well, yielding 

A/’"''(/3,pi,Ti) = ^Pi/''®"“^(T0 + [/(/3+A/3)-/(/3)] 

i 

X - (Pi+Ap*)/(Ti+ATi)l . (22) 


6.2 Applicability 


Our findings pertain specifically to CTI correction employ¬ 


ing the Massey et al. (2014) iterative correction scheme, the 


current nominal procedure for Euclid VIS. Other algorithms 
for the removal of CTI trailing exist that might not be sus¬ 
ceptible in the same way to readout noise. [Prod’homme et al.| 
( |2012[ ), investigating the full-forward approach designed for 
Gaia, did not observe a readout noise floor similar to the 
one we found. The same might hold for including CTI correc¬ 
tion in a forward-modelling shear measurement pipeline (e.g. 
Miller et al.|2013 1. However, the Gaia method has not been 


applied yet to actual observational data, and the |Massey| 
[eni] ( [M4| ) is the most accurate method for the CTI cor¬ 
rection of real data today. 

We remind the reader that our results on the zeropoints 
upon correcting with the correct model (Fig. are depen¬ 
dent on the specifics of the small and faint galaxies we sim¬ 
ulated. Further tests will determine if the large bias in 
persists under more realistic scenarios. 

The narrow tolerances of Apjp = 0.11% and At/t = 
0.17% for the density of the slow traps species might 
look daunting, but fortunately, due to the discernible trails 
caused by these traps it is also the easiest species of which 
to determine the properties. Conversely, the Ap/p = 3% and 
Ar/r = 8% for the fast traps are much larger, but con¬ 
straints on these traps will be harder to achieve from labo¬ 
ratory and in-flight calibration data. Considering the “first 
pixel matched” case, taking into account how trap parame¬ 
ters are determined from CTI trails, relaxes the tolerances 
from ellipticity but tightens the (much broader) tolerances 
from the photometric, for our particular baseline trap mix. 
We notice that, while trap parameters are degenerate and 
Sect. |5.4T^ marks a first attempt to disentangle these param¬ 
eters, each (degenerate) set of parameters can yield a viable 
CTI correction. Characterising the true trap species, how¬ 
ever, is crucial with respect to device physics applications. 

Source profile-dependence of the CTI-induced flux bias 
AFrei will lead to a sample of realistic sources (i.e. with a dis¬ 
tribution of source profiles) showing a range in AFl-ei at any 
given readout noise level. Thus, the uncertainty in Af).ei will 
be larger than the we measured for our broad-winged, 
but homogeneous images in Sect. |5.5.1| More sophisticated 
simulations are necessary to assess the role of the variable 
CTI-induced flux bias in Euclid’s photometric error budget. 


7 CONCLUSIONS AND OUTLOOK 

The goal was to bridge the divide between engineering mea¬ 
surements of CTI, and its degradation of scientific measure¬ 
ments of galaxy shapes. We have developed a very fast algo¬ 
rithm to model CTI in irradiated e2v Technologies CCD273 
devices, reproducing laboratory measurements performed at 
ESTEC. We take a worst-case approach and simulate the 
faintest galaxies to be studied by Euclid, with a broad¬ 
winged exponential profile, at the end of the mission and 
furthest away from the readout register. Our analysis is hin¬ 
dered by the divergent surface brightness moments of the 
Marsaglia-Tin distribution that the ellipticity components 
follow. We alleviate this problem by means of a Taylor ex¬ 
pansion around the mean of the denominator, yielding an 
accuracy of aCa by averaging over 10^ simulated im¬ 

ages. We advocate that Euclid requirements be re-defined 
in a way that avoids ratios of noisy quantities. 

Our detailed study of the trapping process has con¬ 
firmed that not all traps are equally bad for shape measure¬ 
ment ( [Rhodes et al.|[?010[ ): Traps with release timescales of 
a few clocking cycles cause the largest spurious ellipticity, 
while all traps with longer n yield the strongest flux bias. 

The impact of uncertainties in the trap densities pi and 
time scales n on CTI effects can be predicted to a satis¬ 
factory accuracy by a model that is linear in the pi and 
additive in the effects of different trap species. For future 
applications, this will allow us to reduce the simulational 
effort in CTI forecasts, calculating the effect of trap pixels 
from single species data. 

Informed by laboratory data of the irradiated CCD273, 
we have adopted a baseline trap model for Euclid VIS fore¬ 
casts. We corrected images with a trap model M-(-5M offset 
from the model M used for applying CTI. Thus we derived 
tolerance ranges for the uncertainties in the trap parame¬ 
ters, given Euclid requirements, positing that the required 
level of correction will be achieved. We conclude: 

1. In the absence of readout noise, perfect CTI correction 
in terms of ellipticity and flux can be achieved. 

2. Given the nominal rms readout noise of 4.5 electrons, we 
measure = Aei = —1.18 x 10~^ after CTI correction. This 
still exceeds the Euclid requirement of |Aei| < 1.1 x 10“^. 
The requirement may still be met on the actual ensemble 
of galaxies Euclid will measure, since we consider only the 
smallest galaxies of 5/^=1!. Likewise, in S/N = 200 stars, 
we measure a size bias of 1.12 x 10“^, exceeding the require¬ 
ment of I Ai7^/i7true| <4 X 10“"^. 

3. The spurious ellipticity Aei sensitively depends on the 
correct well fill power P, which we need to constrain to an ac¬ 
curacy of A)?toi = (6.31 ± 0.07) X 10“® to meet requirements. 
This assumes calibration by a single, bright charge injec¬ 
tion line. The narrowest tolerance intervals are found for the 
dominant slow trap species in our baseline mix: Aptoi/po = 
(±0.0272 ± 0.0005)%, and Artoi/ro = (±0.0400 ± 0.004)%. 

4. Given the nominal rms readout noise, we measure a flux 
bias Zf = ATrei = (—1.980 ± 0.012)% after CTI correction, 
within the required |Af).ei| <0.7 % for the photometric un¬ 
certainty. More relevant for Euclid will be the uncertainty of 
this bias, which for realistic sources depends on their source 
profile. Further study is necessary here, as well as for the 
impact of CTI on photometric nonlinearity. 
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The final correction will only be as good as on-orbit 
characterisation of physical parameters such as trap loca¬ 
tions, density and release time. The next steps building on 
this study should include: 1.) Researching and testing novel 
algorithms mitigating the effects of read noise as part of 
the CTI correction. 2.) Characterising the effect of realistic 
source profile distributions in terms of the photometric and 
nonlinearity requirements. 3.) Translating the tolerances in 
trap model parameters into recommendations of calibration 
measurements and their analysis, based on modelling the 
characterisation of trap species. 

Plans for Euclid VIS calibration have already been up¬ 
dated to include charge injection at multiple levels such that 
P does not need to be extrapolated from bright charge in¬ 
jection lines to faint galaxies. We will continue to liaise be¬ 
tween engineers and scientists to determine how accurately 
it will be necessary to measure these physical parameters. 
The VIS readout electronics will be capable of several new 
in-orbit calibration modes such as trap pumping ( [Murray] 
et al.|[Ml2| that are not possible with HST, and our calcu¬ 


lations will advise what will be required, and how frequently 
they need to be performed, in order to adequately charac¬ 
terise the instrument for scientific success. 
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pixels 



pixels 

Figure Al. CCD273 EPER trails in the serial {upper plot) and 
parallel {lower plot) directions. Shown here are the G quadrant 
trails at an input signal of ~ 1000 electrons. Solid lines within 
the light and dark grey shaded areas denote the average and its 
uncertainty of the profile before and after correction for electronic 
effects. The best-fit model to the corrected trail is shown as a 
long-dashed line. For the purpose of illustration, the baseline trap 
model is shown in both plots as a short-dashed line. Building on 
the serial EPER model, the baseline model includes fast traps 
that are seen in quadrant F. 


APPENDIX A: INFORMING THE BASELINE 
MODEL WITH LABORATORY DATA 


Al EPER/FPR data with irradiated CCD 


In this Appendix, we define a baseline CTI model for Eu¬ 
clid VIS. Our model is based upon laboratory tests of an 
irradiated e2v Technologies back-illuminated Euclid pro¬ 
totype CCD273, analysed at ESA/ESTEC jProd’homme 
|et al.|[2Q14[ ). The device was irradiated at ambient room 
temperature using 10.4 MeV protons, degraded from a 
38.5 MeV primary proton beam at the Kernfysisch Ver- 
sneller Instituut, Groningen, in April 2013. Two different 


shielding masks were used (Prod’homme et al. 20141 re¬ 


sulting in the four quadrants of the CCD, called E, F, G, 
and H, and corresponding to the four output nodes, re¬ 
ceiving different radiation doses. Each a half of two quad¬ 
rants, called G and H, received a 10 MeV equivalent flu- 
ence of 4.8 x 10® protons/cm“^, representative of the pre¬ 
dicted end-of-life (eol) proton fluence for Euclid. Half of the 



Signal level [electrons] 



Signal level [electrons] 


Figure A2. The well fill power /3 measured from the integrated 
EPER CTI as a function of input signal. The upper panel shows 
the results from the serial EPER measurements, for which CTI is 
present in the F and G quadrants and can be corrected using the 
E and H quadrants. The lower panel shows the results from the 
parallel EPER measurements, for which CTI is present in the F, 
G, and H quadrants and can be corrected using the E quadrant. 
Open symbols denote the raw measurements, filled symbols the 
calibrated measurements from which the fits for 0 are derived. 


F quadrant was irradiated with a 10 MeV equivalent flu¬ 
ence of 2.4 X 10® protons/cm“®, the eol/2 fluence. Neither 
the E quadrant, the serial register of the H quadrant, nor 


the readout nodes were irradiated {[Verhoeve et al. 2014 
Prod’homme et al.|2014|. 


At the ESA Payload Technology Validation section 
GCD test bench located at ESTEG [ Verhoeve et al.||2014 |, 
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the irradiated CCD273 was characterised at the Euclid VIS 
nominal conditions of 153 K temperature and a 70 kHz read¬ 
out frequency. While a serial clocking scheme with the same 
width for each register phase at each step was used, minimis¬ 
ing serial CTI, the nominal line/parallel transfer duration of 
0 .11ms was not optimised. 

As part of the characterisation, a suite of extended pixel 
edge response (EPER) and first pixel response (FPR) ex¬ 
periments were performed, at different flatfield signal levels. 
For the purpose of deriving a fiducial baseline model of the 
charge traps present in the CCD273, we focus on the parallel 
and serial EPER data. To study the serial EPER (sEPER) 
CTI, a flatfield image is taken, then the half opposite to the 
readout direction is dumped; then the frame is read out. 
This yields a flatfield with a sharp trailing edge in flatheld 
signal. Electrons captured from flatheld pixels are being re¬ 
leased into signal-less pixels, resulting in a CTI trail. Our 
parallel EPER (pEPER) tests make use of the parallel over¬ 
scan region, providing a similar signal edge. 

Each measurement was performed repeatedly, in order 
to gather statistics: 45 times for the sEPER data at low sig¬ 
nal, and 20 times for the pEPER data. Raw trail prohles 
are extracted from the hrst 200 pixels following the signal 
edge, taking the arithmetic mean over the independent lines 
perpendicular to the direction (serial or parallel) of interest. 
The same is done in the overscan region, unexposed to light, 
and the pixel-by-pixel median of this reference is subtracted 
as background from the raw trails. In the same way as the 
reference, the median hatheld signal is determined, and also 
corrected for the overscan reference. Finally, the trail (hat- 
held signal) at zero hatheld exposure time is subtracted from 
the trails (hatheld signals) at exposure times >0. 

Figure f^T] shows the resulting “uncalibrated” trail pro¬ 
hles for the sEPER (upper panel) and pEPER (lower panel) 
measurements in the G quadrant (eol radiation dose), at a 
hatheld exposure time corresponding to an average of 1018 
signal electrons per pixel. These are the upper solid lines 
with light grey shading denoting the propagated standard 
errors from the repeated experiments. Effects in the readout 
electronics mimic CTI. We correct for the electronic effect by 
subtracting the average trail in the unirradiated quadrants 
(E for pEPER, and E and H for sEPER). The resulting “cal¬ 
ibrated” trail prohles and their uncertainties are presented 
as the lower solid lines and dark grey shadings in Fig. |A1[ 
The calibration makes a small correction to the sEPER trail 
which is dominated by slow traps, yielding a signihcant sig¬ 
nal out to ~60 pixels. On the contrary, the electronic ehect 
accounts for 1/3 of the uncalibrated pEPER trail even in 
the hrst pixel, and for all of it beyond the tenth. Thus the 
S/N in the calibrated trail is much lower. 


A2 The well fill power /3 

In a volume-driven CTI model, the cloud of photoelectrons 
in any given pixel is assumed to hll up a height within 
the silicon that increases as electrons are added (Eq. 121. 
The growth of the clou d volume is governed by the term 
in Eh- (12l, with the full-well depth w = 84700 
limiting the maximum number of electrons in a pixel. There 
is no supplementary buried channel in the CCD273, which 


Table Al. The same as Table[^ but for the best-fit models shown 
in Fig. The baseline well fill power is /3o = 0.58. 


best-fit sEPER model 
Trap density pi [px“^] 
Release timescale ti [px] 

i = l 

0.01 

0.8 

i = 2 

0.03 

3.5 

i = 3 

0.90 

20.0 

best-fit pEPER model 

i = l 

i = 2 

i = 3 

Trap density pi [px“^] 

0.13 

0.25 

— 

Release timescale [px] 

1.25 

4.4 

— 


for HST/ACS leads to the first ~ 100 electrons effectively 
occupying zero volume ( [Massey et al.|[^10[ ). 

We use measurements of the integrated EPER as a func¬ 
tion of input signal to constrain the well fill power /3 of the 
trapping model. Our simulated galaxies are faint; so we re¬ 
strict ourselves to the four lowest signal levels measured in 
the laboratory, with up to ~ 1000 electrons. The input sig¬ 
nal is measured as the average count difference between the 
flatfield and overscan regions, corrected for the CCD gain. 

Figure |A2| shows the CTI trails from Fig. |A1[ inte¬ 
grated over the first 12 pixels. We checked that integrat¬ 
ing over up to the full overscan region length of 200 pix¬ 
els does not change the results drastically. In the sEPER 
data (upper panel of Fig. A21, the unirradiated quadrants E 
and H (open squares and diamonds) exhibit very small trail 
integrals (caused by the readout electronics); one order of 
magnitude smaller than in the irradiated quadrants F and 
G (open circle and triangle). Hence, calibrating out the in¬ 
strumental effect by subtracting the arithmetic average from 
the E and H quadrants yields only a small correction to the 
F and G trail integrals. To these calibrated sEPER mea¬ 
surements (filled circle and triangle), we fit linear relations 
in log-log-space using the IDL fitexy routine and measure 
/Sp.cai = 0.49 ± 0.04 and ^dc.cai = 0.58 ± 0.03. 

We repeat this procedure for the pEPER measurements 
where the unirradiated E quadrant shows a similar EPER 
integral than the irradiated F, G, and H quadrants (lower 
panel of Fig. A2|. Thus, the pEPER and sEPER integrals 
may yield similar values as a function of signal, but for 
pEPER the low S/N<^ 1 causes large uncertainties. Conse¬ 
quently, P is not well constrained, with /lp,cai = 0.66 ± 0.53, 
/3G,cai = 0.61 ± 0.36, and /3H,cai = 0.61 ± 0.89, but they agree 
with the sEPER results. 

In conclusion, we adopt a baseline well-fill power of Po = 
0.58 for our further tests, based on the precise sEPER result 
for the full radiation dose. 


A3 Fhom trail profiles to trap parameters 


To constrain the trap release time-scales Tt and trap densi¬ 
ties pi , we make use of the two signal levels of ~ 360 electrons 
and ~ 1000 electrons that bracket the number of electrons we 
expect to be found in a typical faint Euclid galaxy. These are 
the two highest data points in Fig. |A2| We compare the av¬ 
erage, measured, calibrated trails from the irradiated quad¬ 


rants (examples for the G quadrant are presented in Fig. Al I 
and compare them to the output a one-dimensional version 
of our [Massey et al.j ( |2014[ ) clocking routine produces given 
trap densities pi and release timescales ri, and under circum- 
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Figure A3. Trail profiles corresponding to our trap models. Upper left panel: Trails for models in which the densities pi are modified 
by a factor 0.9 (red) or l.f (blue). This is applied to: the fast traps pi (dashed lines); the medium traps p 2 (long-dashed lines); the slow 
traps p 3 (dot-dashed lines); pi and p 2 (triple dot-dashed lines); all traps (dotted lines). Upper right panel: The same, but relative to 
the baseline model. Lower left panel: Trails for trap models in which the release time scales Ti are modified by a factor 0.7 (red) or 1.4 
(blue). The line code is the same as above. Solid lines show models in which not only all the r; are changed, but also the pi adjusted 
such that the first pixel matches the baseline trail profile (Sect. [^.4.4[). Lower right panel: The same, but relative to the baseline model. 


stances close to the laboratory data (i.e. a 200 pixel over¬ 
scan region following a 2048 pixel flatfield column of 1018 
signal electrons). The fitting is performed using the MPFIT 
implementation of the Levenberg-Marquardt algorithm for 
nonlinear regression ( Markwardt|2009 More|[l978 1. 


Fitting a sum of exponentials is remarkably sensitive to 
noise in the data because the parameters {n and pi) we are 
probing are intrinsically correlated. We assess the robustness 
of our results by repeating the fit not only for the two (three) 
irradiated sEPER (pEPER) quadrants at two signal levels, 
but for a wide range of trail lengths (60 < K < 150) we 
consider, and with and without adding a constant term. 


There are several possible trap species as defined by 
their n that show up in our data set. We rule out those of 
very low densities and consider the frequent “species” whose 
time-scales are within a few percent of each other as one. 
Still, this leaves us with more than one contesting family of 
trap species that yet give similar trails in some of the quad¬ 
rant/signal combinations. Because, at this stage, our goal is 
to derive a plausible baseline model rather than pinpointing 
the correct trap species, we filter for the most common Ti and 
give precedence to the higher-S/N data (sEPER, end-of-life 
dose, 1000 signal electrons). The resulting best-fit models are 
shown in Table [AT] and Fig. |A1| The actual baseline model 
(Table short-dashed line in Fig. |A1[ | includes additional 
fast traps seen in the lower-S/N data. We raise the density 


from 0.94 traps per pixels to a mnemonic total of 1 trap per 
pixel at end-of-life dose. More refined methods will be used 
to determine the trap species in a more detailed analysis of 
irradiated CCD273 data. 

Because only 464 pixels of the serial register in the test 
device were irradiated, the effective density of charge traps 
an electron clocked through it experiences is smaller by a 
factor of 464/2051 than the actual trap density correspond¬ 
ing to the end-of-life radiation dose that was applied. We 
correct for this by quoting results that have been scaled up 
by a factor of 2051/(464 x 0.94)«4.155. 


A4 Example CTI trails 

Figure [A3| shows, for the largest deviations from the baseline 
trap model we consider, their effect on the shape of the CTI 
trails. Using our CTI model, we simulated the trail caused 
by a single pixel containing a signal of ~ 1000 electrons, 
comparable to a hot pixel in actual CCD data. 
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